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Abstract 

We consider the period-doubling and intermittency transitions in 
iterated nonlinear one-dimensional maps to corroborate unambigu- 
ously the validity of Tsallis' non-extensive statistics at these critical 
points. We study the map x n+ \ = x n + u \x n \ z , z > 1, as it describes 
generically the neighborhood of all of these transitions. The exact 
renormalization group (RG) fixed-point map and perturbation static 
expressions match the corresponding expressions for the dynamics of 
iterates. The time evolution is universal in the RG sense and the non- 
extensive entropy Sq associated to the fixed-point map is maximum 
with respect to that of the other maps in its basin of attraction. The 
degree of non-extensivity - the index Q in Sq - and the degree of non- 
linearity z are equivalent and the generalized Lyapunov exponent X q , 
q = 2 — Q~ l , is the leading map expansion coefficient u. The corre- 
sponding deterministic diffusion problem is similarly interpreted. We 
discuss our results. 

1 Introduction 

The proposal, in 1988, by Tsallis [1] of a non-extensive generalization of sta- 
tistical mechanics unleashed an unprecedented discussion [2] on the founda- 
tions of this most basic branch of physics. This is an exceptional event, since 
its inception, with Boltzmann and Gibbs (BG), in the long and heretofore 
unchallenged history of statistical mechanics. Many recent studies [2] have 
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offered experimental and numerical evidences indicating limiting conditions 
for the domain of validity of the canonical BG theory. But, so far there have 
been no firm developments to tilt the balance towards universal acceptance 
of these new ideas, nor a definite understanding of the deep-lying physical 
reasons (believed in some cases to be a breakdown in the chain of increasing 
randomness from non-ergodicity to completely developed chaoticity) for the 
failure of the BG statistics and the competence of the non-extensive gener- 
alization. However, here we present exact analytical results for a class of 
critical points in non-linear maps from which a rigorous corroboration and 
understanding of the non-extensive theory can be drawn [3] - [7]. 

The renowned period-doubling and intermittency routes to chaos are now 
classical examples of how the use as starting points of simple non-linear dis- 
crete maps have often led to important developments in the theory of non- 
linear dynamical systems [8]. The discovery of the universal properties asso- 
ciated to these mechanisms, comparable with those of conventional critical 
phenomena in statistical physics, triggered, over two decades ago, an explo- 
sion of activity in the field and now both routes are well understood and 
verified experimentally. Here we might add yet other universal aspects to 
their properties (explained in more detail below) the non-extensivity marked 
by the power-law growth or decay of the sensitivity to initial conditions, and 
the entropy extrema associated to each of the fixed points that define the 
RG universality classes. 

The well-known period-doubling and intermittency transitions are based 
on the pitchfork and the tangent bifurcations. In particular, the one-dimensional 
logistic map and its generalization to non-linearity of order 3 > 1, f^x) = 
1 — /i |x| 3 , — 1 < x < 1, exhibit these two different types of bifurcations an in- 
finite number of times as the control parameter reaches critical values. At 
these critical points universal scaling laws and indexes hold, independently 
of the details of the map, and are therefore shared by all maps of the same 
order 3. As with other statistical mechanical systems an explanation of uni- 
versality is provided via the renormalization group (RG) method, which has 
been successfully applied to this type of iterated map. The RG doubling 
transformation, consisting of functional composition and rescaling, was first 
devised to study the cascade of period-doubling transitions and its accumula- 
tion point, and was later applied to the intermittency transition. As pointed 
out [8] the latter case is one of the rare examples where the RG equations 
can be solved exactly. Here we indicate that this analytic solution holds 
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generally not only for tangent bifurcations (as known time ago [8]) but also 
for pitchfork bifurcations (as noticed recently [4]). We also point out that 
the fixed-point map can be determined as an extremum of the non-extensive 
entropy, thus providing evidence for a possible connection between extremal 
properties of entropy expressions and those for the fixed points obtained from 
the usual RG recursive method [9]. 

As required, at each of the map critical points the Lyapunov Ai exponent 
vanishes, and the sensitivity to initial conditions £ t for large iteration time t 
ceases to obey exponential behavior, exhibiting instead power-law behavior 
[10]. As the means to describe the dynamics at such critical points, and 
based on the non-extensive entropy of Tsallis Sq [2] (see Eq. (16) below), 
the g-exponential expression 

& = exp 9 (V) = [!-(?- 1) V]" 1/(ff_1) (!) 

containing a g-generalized Lyapunov exponent \ q was proposed [11], together 
with generalizations for the rate of entropy production (referred here as the 
Kolmogorov-Sinai (KS) entropy [8]) Kq (see Eq. (15) below) and for the 
Pesin identity X q = Kq [11]. The standard expressions are recovered when 
q, Q — > 1. The use of the two indexes q and Q in our application will become 
clear below. Furthermore, recent studies [11]- [14], focusing mainly at the 
onset of chaos (the accumulation point of the period-doubling bifurcations) 
of the 3-logistic map have revealed a series of connections between the Tsallis 
entropic index q and the non-linearity 3, the fractal dimension df of the 
chaotic attractor, and the end points of its multifractal singularity spectrum 
f(a). Here we identify q and \ q in terms of the tangent and pitchfork critical 
map parameters. 

The properties of dynamical systems at critical points can also be probed 
by considering the deterministic diffusive processes that take place in non- 
linear maps with discrete translational symmetry. These maps consist of a 
sequence of cells containing each elements of a basic nonlinear map. One such 
construction [15] makes use of the tangent bifurcation, and, interestingly, 
anomalous diffusional behavior has been observed close to the intermittency 
transition [15], [16]. Continuous-time random- walk theory concepts have 
been applied [15], [16] and, characteristically, the associated waiting-time 
distribution function ip(t) has been shown to have non-exponential power- 
law decay. The validity of the non-extensive approach to criticality has been 
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recently checked by examination of the anomalous diffusion at the intermit- 
tency transition [16], [17]. The key results in these studies center on the 
analytical form predicted for ip{t) which bears a close relationship to the 
power-law Eq. (1). Here we point out that the analytical forms for this 
and related distributions actually correspond to those obtained by the re- 
peated application of the RG transformation. As mentioned, if the fixed 
point solution can also be determined via an entropy optimization method 
the expressions for these distributions are indicative of the same entropy 
extremum. 

2 Outline of results 

The main points in the following analysis are: 

1) The connection between the RG fixed-point map expression for the 
iterate x t associated to trajectories x t with xo — > at t = and the expression 
in Eq. (1) for the sensitivity of initial conditions £ t . This leads to simple 
identifications of the non-extensive parameters as q = 2 — z^ 1 and X q ~ a. 
We corroborate the earlier known results q = 3/2 and q = 5/3 for the tangent 
and pitchfork bifurcations for the logistic map 3 = 2, respectively [11]. 

2) The connection between the RG fixed-point map recursion relation for 
the iterate x t and the expression obtained from Sq for the temporal evolution 
of the number of cells W t occupied by an ensemble of initial conditions spread 
over a finite-length interval [0, Xo]. (See the definition of Sq in Eq. (16), and 
the explanatory comments at the end of this section about the relationship 
q = 2 — Q~ x between the indexes q and Q). This leads to the identifications 
Q = z and K q ~ a, confirming q = 2 — z^ 1 and giving support to the validity 
of the generalized Pesin identity X q = Kq at the critical transitions. See Ref. 
[7] for a recent proof of the validity of this generalized identity at the edge 
of chaos of the logistic map 3 = 2. 

3) The general property of the non-extensive entropy function which de- 
creases monotonically as the RG transformation flows away from the fixed 
point where it attains its maximum value. 

4) The maximum entropy property of the distribution ^(r), ^(r) = 
Jq dt ip{t), of cell residence time intervals [0, r] in the critical diffusive pro- 
cesses. 

For clarity of presentation it is convenient to advance the following prop- 
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erty of the g-exponential function. For this function the ordinary derivative 
dexp(x)/dx = exp(x) and inverse [exp(a;)]~ 1 = exp(— x) properties become, 
respectively, dexp q (x)/dx = [exp q (x)] q and [expire)] -1 = exp 2 _ (J (— x), and 
these combine into the identity 

Q 

(2) 

where q = 2 — Q^ 1 . Thus, in a given problem where the g-exponential 
function applies there appear related values of the g-index, that become the 
same (and the above identity becomes trivial) as q — > 1. 

3 Critical dynamics from fixed-point map 

So, we recall [8], [19] the solution of Hu and Rudnick to the Feigenbaum RG 
recursion relation obtained for the case of the tangent bifurcation. For the 
transition to periodicity of order n consider the n-th composition of the 
original map / in the neighborhood of one of the n points tangent to the line 
with unit slope and shift the origin of coordinates to that point. Then, one 
obtains 

f( n )(x) = x + u\x\ z + o(\x\ z ), (3) 

where u > is the leading expansion coefficient (the right- hand- side in Eq. 
(3) may depend on n only through the omitted terms). The RG fixed-point 
map x' = f*(x) was found to be 

x' = xexp^ux*- 1 ) = x[l -(z- ljua;*- 1 ]- 1 /**- 1 ), (4) 

as it satisfies /*(/*(#)) = a~ l f*(ax) with a = 2 1 ^ z ^ and has a power- 
series expansion in x that coincides with Eq. (3) in the two lowest-order 
terms. (We have used x z ~ x = |a;| z_1 sgn(rr)). Eq. (4) can be rewritten as 
a;'-^- 1 ) — x~( z ~ lS) — (z — l)u, a translation of the power law y = x~^ z_1 \ 
the property that led Hu and Rudnick in their derivation. The effect of a 
perturbation of the form x~ v transforms this into 

x'- (z -^ + e'x'- p = x^ z - 1] + ex- p - (z - l)u, (5) 

or, to lowest order in e' = e, 

x' = xexp z (-ux 2 " 1 ) [exp z (-ua; z_1 )] 



exp (x) = 



exp (|) 
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z—p 




z-l\P~P 



(6) 
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The recurrence relation in Eq. (5) remains invariant under iteration and 
rescaling by a when e' = e is multiplied by the factor 7 = 2( p ~ z+1 ^( z ~ 1 \ As 
it is known [18], [19], the leading relevant perturbation about this critical 
point is that which corresponds to a shift from tangency, for this case the 
lowest order term in x in Eq. (6) must be of order one, and so we fix p = z 
for the perturbation exponent (in the unormalized eigenf unction). 

Fittingly, it can be readily concluded that the RG recursion relation in 

the neighborhood of the pitchfork bifurcations df^ 1 \x)/dx = —1 con- 

x=0 

forms with the above scheme and that its solution is also applicable to these 
transitions provided the sign of u is changed for x > 0. It is interesting to 
point out that at the time the RG method was first applied to the tangent 
bifurcation it was thought that the simple closed form obtained for the fixed 
point map, the g-exponential in Eq. (4), would most likely prove not to be 
fruitful for the study of period doubling [8], [19]. Not only is this solution ap- 
plicable to each period doubling transition with the minor modification of the 
indicated sign change, but also, as it has been recently demonstrated [6], the 
time evolution at their accumulation point, the edge of chaos, is expressable 
in terms of similar g-exponentials. Below we comment on the caveats that 
appear in relating z above to the 3 in the initial map e.g. = 1 — fi \x\ 3 . 

When the number of iterations (after an initial Xq) is large and the con- 
tinuous time t approximation is taken ( X* — Xf) Xq — X and u — at in Eq. 
(4)) the fixed-point map can be written as 



with properties dx t /dt = ax z and dx t /dxo = (x t /x ) z . Thus, the fixed-point 
and perturbation relations can be used to describe the time evolution of 
intervals containing iterates, an initial interval [0, xo < x + ] of length Axq 
evolves into [0, x t < x + ] with length Ax t where x t is the image of x at time 
t. Accordingly, Eq. (5) with (p = z) becomes 



We should like to emphasize that the expression for time evolution Eq. (7) 
has the same universal form as that for the fixed-point map Eq. (4). In the 



x t = x [l-(z- 1)qXq 




(7) 



Ax;^ + e t Axt z = Ax (z - 1] + e Ax z - (z - l)at. 



(8) 
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neighborhood of the points of tangency f^ n \0) the iterates follow monotonic 
paths with time shape set by the map shape itself, i.e. long time dynamics 
follow the static solution of Hu and Rudnick. It should also be kept in 
mind that Eq. (8) applies only to small perturbations as this corresponds 
to the perturbation recursion relation obtained to linear order around the 
RG fixed point. We also comment that our use here of the continuous time 
approximation (valid for large iteration time t) in deriving Eq. (4) is only 
for expediency reasons and has no effect in the generality of our results as 
this equation is equally obtained from the scaling property 

x m = r™(x)=m- 1 '*- 1 r(m 1 '*- 1 x), m = l,2..., (9) 
that in turn is obtained from f*(f*(x)) = a^ 1 f*(ax) with a = 2 1 ^ z ~ 1 [5]. 

4 Sensitivity to initial conditions 

We have reached a point at which a straightforward proof of the validity of 
Tsallis' Eq. (1) can be obtained. This is done by direct evaluation of the 
sensitivity to initial conditions from its definition £ t = lim^ xo _, (Ax t / Axq). 
From Eq. (7) and dx t /dx = (x t /x Q ) z one obtains 

&x ) = [1 - (z - VaxS-H]-'"'- 1 ), (10) 

and comparison with Eq. (1) yields q = 2 — z~ x and \ q (x ) = zax^lW}. 

We can also confirm the values of q obtained numerically by Tsallis et 
al [11] for both the tangent and pitchfork bifurcations when 3 = 2. The 
composition fffi for the 3-logistic map f^x) = 1 — fj, \x\ 3 , 3 > 1, at the 
tangent bifurcations is fj?\x) = x + ux 2 + o(x 2 ), u > 0, and from z = 2 
one has q = 3/2. For the pitchfork bifurcations one has instead fjT\x) = 
x + ux 3 + o(x 3 ), because d 2 f^ ^ jdx 2 = at these transitions and u < is 
now the coefficient associated to d jdx < 0. In this case we have z = 3 
in q = 2 — z^ 1 and one obtains q = 5/3. Notably these specific results for 
the index q are valid for all 3 > 1 and therefore define the existence of only 
two universality classes for unimodal maps, one for the tangent and the other 
one for the pitchfork bifurcations [5]. 

Further, if one uses the fact [20] that the invariant distribution of 
in Eq. (3) is of the form p(x) ~ x~^ z_1 \ the average \ q of \ q (x ) over x 
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yields the constant X q ~ a. In relation to this we observe that this average 
is compatible with 

the g-extension of the customary [8] expression for Ai as the average of 
In \df*(x)/dx\ over p(x). By taking the g-logarithm (where ln q y = (y 1 ~ q — 
1)/(1 — q) is the inverse of exp (?/)) of both sides of 

Ux ) = [l-(q-l)\ q (x )t]- 1 ^- 1 \ (12) 
and by recalling the definition of £t(^o) one obtains 

X q (x ) = r l In, 

where ln g \ df*^ (x ) / dx \ ~ t as X q (x ) does not depend on t. (Notice that 
the limit t — > oo of Eq. (13) is the g-extension of the regular definition of 
Ai(x ) [8])- Next, use of the fixed-point relation /* ( % ) = i 1 / 1- */*^ 1 '* -1 ^) 
andt = (x/^o) 2 " 1 into (13) leads to X q (x ) = (x/x )~( z ~V \n g \df*(x)/dx\ and 
Eq. (11). 



df*(x) 



dx 



(11) 



df^{x Q ) 



dxo 



(13) 



5 Rate of entropy change 

Even though our basic results have been derived, it is useful to extend our 
analysis and reveal more details about the dynamics at these critical points. 
Our following arguments focus on x t rather than dx t /dxo and imply Q = z, 
thus confirming q = 2 — z~ x . They also indicate that the rate of entropy 
change is Kq ~ a, suggesting the validity of the g-generalization of the Pesin 
\ q = Kq at these critical states. Therefore, let's consider a cell partition 
(N — * oo cells of length I — > 0) in the 'entire' phase space [0, £+]. Consider 
also the time evolution of occupied cells W t by an ensemble of iterates with 
an extended initial occupation W that spans [0, x < x + ]. We write the cell 
distribution of iterates pi as 

1 _|_ r W 

^ = ^A5> = °> ( 14 ) 



8 



as this is a convenient form when considering normalized perturbations to the 
uniform distribution. It is a general expression when no further conditions 
are imposed on the 5{ but below we require small S^. When Eq. (14) is 
introduced into the KS g-entropy difference 

K Q t = S Q (t)-S Q (0), t large, (15) 

where 

l- y w P Q 

s Q = ^r\ Vi (i6) 



Q-l 



is the Tsallis entropy, we obtain 



W t Wo 

W t ~ Q £(1 + 5 ht f = W Q £(1 + 5 ij0 ) Q - (Q - l)K Q t. (17) 

i i 

To lowest order in this reduces to 

w7 (g_1) + \q{q - mw t - Q = 

W - {Q ~ l) + \q(Q - 1)5 2 W - Q - (Q - l)K Q t, (18) 
with 5 2 = Ef si 

Remarkably, Eqs. (18) and (8) can be seen to be equivalent if Q = z and 
if W t fi ~ Ax tj o, Kq ~ a and |Q(Q — l)£ t 2 ~ £ t,o- The entropic index Q 
is simply given by the non-linearity z while Kq is the expansion coefficient 
a (and through the proviso that X q = Kq this would imply X q ~ a). Our 
comparison between Eqs. (18) and (8) is consistent in that each is the lowest- 
order expression valid for small perturbations, the first to the uniform cell 
distribution of iterates and the second to the fixed-point map solution. We 
can draw the following conclusion from the association of Eqs. (18) and 
(8): The entropy Sq for the fixed-point map is maximum, as Sq — Sq = 
Q5 2 /2W q > 0, and the distribution becomes uniform, pi = W~ x . The RG 
flow is away from the fixed-point since 7 > 1 and the entropy decreases as 
the RG transformation is applied. The proposed equivalence between Eqs. 
(18) and (8) is justified by making the following connection with Eq. (1). 
For a uniform distribution Wtfi ~ Axtfi, as the ratio of the number of cells 
W t occupied at time t to the distance Ax t is the length /, and in the limit 
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N — > oo and I — > 0, so that Axo — > 0, with say Axq/1 = b, one obtains 
W t = b£ t - This explicit link between the occupation numbers W t and £ t 
is in line with the all-important connection between the loss of information 
embodied by the KS entropy and the sensitivity to initial conditions in non- 
critical regimes [8], and with the parallel connection suggested by Tsallis et 
al [11], [12] for the analogous quantities at critical points. 

It is worth the mention that the use of g-exponential (and the inverse q- 
logarithm) functions in the derivation of the properties of a given system will 
in general introduce pairs of conjugate indexes such as q = 2 — Q^ 1 above. 
As pointed out, taking derivatives and inverse operations lead to expressions 
involving one or the other of these simply related values of the entropic 
index (that become trivially the same q = Q = 1 for the usual exponential 
and logarithmic functions). Consequently, while some theoretical features 
are expressed through q some others appear expressed via Q, but in all cases 
their occurrence and meaning is apparent. 

6 Deterministic diffusion 

We comment now on the deterministic diffusion process that takes place for 
the variable x t in the map x t+1 = f(x t ), where f(m+x) = m+f(x), m integer 
and f(—x) = —f(x). Geisel and Thomae [15] found that when the unit-cell 
map /(0 < x < 1/2) is chosen to have the (intermittency transition) form 
x t+ i — x t + ax\ within [0,x + < 1/2] anomalous diffusion, of the dispersive 
type, is generated when z > 2 since then the mean-square displacement 
< Ax 2 (t) >~ t 1 /^ -1 ). This behavior follows [15] after the distribution of cell 
residence times ip{t) is found to be 

ip(t) = ax^l + (z- l)ax z -HY z/{z -^. (19) 

Grigolini and co-workers [16] recognized the relationship between the ex- 
pressions for ip(t) and They also analyzed, and verified, the validity of 
the non-extensive scheme involving the g-generalized KS entropy its equality 
with the g-generalized Lyapunov exponent [16], [17]. Here we can recap these 
results as follows: The properties of the diffusion process driven by the in- 
termittency transition are universal since they stem from the RG fixed-point 
map expression. Namely, ip{t) is, besides normalization, the derivative of the 
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inverse fixed-point map dxo/dxt with and is therefore inversely pro- 

portional to £t = dxt/dxQ. The distribution ^(r) = J T rft ^(i) ~< Ax 2 (t) > 
is in turn proportional to the inverse of the fixed-point map x (x T ). Also, 

*(r) = A[l + (z - ^a^V]- 1 /^- 1 ), (20) 

has the form of the distribution obtained by optimization of the non-extensive 
entropy with index Q = z for fixed ax z + x r. The related entropy index 
q = 2 — Q^ 1 = 2 — applies to the residence time distribution ip(t). 

7 Summary and discussion 

Thus, for the unimodal maps analyzed we have indicated that the exact static 
solution of the RG equations for the tangent and pitchfork bifurcations also 
describe the dynamics of iterates. Further, these fixed-point expressions are 
the same g-exponential expressions describing the temporal evolution of en- 
sembles of iterates prescribed by the non-extensive formalism for non-linear 
maps at criticality. Therefore the RG fixed-point map for the universality 
classes determined by z can be constructed via Tsallis entropic arguments. 
The equivalence rests on the derivation of Eq. (1) exclusively from RG pro- 
cedures. The study of the pitchfork and period-doubling transitions has been 
expanded via detailed derivation of their g-generalized Lyapunov exponents 
X q and interpretation of the different types of sensitivity £ t [5]. Likewise, 
the properties of the intricate trajectories at the edge of chaos of the logistic 
map 3 = 2 have been analytically obtained leading too to the determination 
of X q and interpretation of the dynamics at the strange attractor [6]. Fur- 
ther, the validity of the g-generalized Pesin identity at this important state 
has been recently rigorously proved [7]. More generally, we addressed here 
the properties of the fixed-point map with the purpose of finding additional 
evidence [9] for an existing relationship between the variational properties 
of entropy expressions and the RG approach when applied to systems with 
scale invariance properties. We found that the non-extensive entropy func- 
tion decreases monotonically as the RG transformation flows away from the 
fixed point where it attains its maximum value. 

The ergodic hypothesis lies at the foundation of statistical mechanics im- 
plying that trajectories in phase space cover uniformly the entire pertinent 
regions. But is this hypothesis always correct? Already many years ago the 
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answer to this question has been probed by the study of simple dynamical 
systems with only a few degrees of freedom [8]. Besides their uncomplicated 
description these systems display extremely convoluted motion in phase space 
that is neither regular nor simply ergodic, and the mechanism by which er- 
godicity emerges in these and more complex deterministic systems has been 
effectively explored by studying the sensitivity to initial conditions and the 
associated Lyapunov exponents [8]. The distinction between periodic and 
chaotic motion is signaled, respectively, by the long-time exponential ap- 
proach or departure of trajectories with close initial positions. Here we have 
analyzed the borderline critical states in one-dimensional non-linear maps at 
which the exponential sensitivity law stops working and find that the uni- 
versal dynamical behavior at these states precisely follow the predictions of 
the generalized non-extensive theory. 

The currently developed non-extensive generalization [1], [2] of BG sta- 
tistical mechanics entails a rare examination of the domain of validity of this 
most fundamental base of physics. The suggested physical circumstances 
for which BG statistics fails to be applicable are thought to be associated to 
states that lack the full degree of chaotic irregular dynamics that probes phase 
space thoroughly and that is necessary for true equilibrium. Such anomalous 
states are signalled by the vanishing of the largest Lyapunov exponent and 
exhibit non-ergodicity or unusual power-law mixing. One particular case in 
point is that of critical states in nonlinear unimodal maps. Here impeded or 
incomplete mixing in phase space (the interval — 1 < x < 1) arises from the 
special 'tangency' shape of the map at the pitchfork and tangent transitions 
that produces the monotonic trajectories given by Eq. (7). This has the 
effect of confining or expelling trajectories causing anomalous phase-space 
sampling, in contrast to the thorough coverage in generic states with Ai > 0. 
By construction the dynamics at the intermittency transitions, describe a 
purely nonextensive regime, since the map studied here, Eq. (3), does not 
consider access of trajectories to an adjacent or neighboring chaotic region, 
as in the setting of Refs. [16], [17] or as in trajectories in conservative maps 
with weakly developed chaotic regions [21]. Hence there is no reappearance 
of trajectories from chaotic regions that would cause the relaxation from 
the nonextensive regime with vanishing ordinary Lyapunov exponent to an 
extensive regime with a positive one at some crossover iteration time r. 

The links we have exhibited through exact analytical results between the 
various properties of one-dimensional non-linear maps at critical conditions 
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provide a clear-cut corroboration, with a universal attribute, in the RG sense, 
of the validity of non-extensive statistics at such states. 
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